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Abstract 


We present an analytical result for the supercurrent across a 
superconductor/quantum-dot/superconductor junction. By converting the 
current integration into a special contour integral, we can express the cur- 
rent as a sum of the residues of poles. These poles are real and give a natural 
definition of the Andreev bound states. We also use the exact result to explain 
some features of the supercurrent transport behavior. 

PACS numbers: 74.50.+r, 73.40.Gk, 73.20.-b, 73.63.Kv. 


Introduction. The mesoscopic superconductor/quantum-dot /superconductor (S-QD-S) 
junction is a typical structure to study the phase coherent current transport in mesoscopic 
hybrid systems [II]. When using Keldysh Green function technique |P] to evaluate the current 
across the S-QD-S junction, the usual method demands a numerical integration and faces 
the problem of explaining some features of the current transport behavior [B]. In this paper, 
we employ a new method to analytically compute the current across the S-QD-S junction in 
the absence of voltage bias, and use the analytical results to explain some features of current 
transport properties. 

Model and Hamiltonian. In the system under consideration, a quantum dot defined in a 
2-dimensional electron gas (2DEG) is coupled to two BCS superconducting leads. We model 
this system by the Hamiltonian 


H=H,+Ha+Hprt+Hr, (1) 


where 


Hr, = 5 eval ako + Y.[Arajyal yy + Aj akart] (2) 
k,o k 

Hy = oldd, + did) , (3) 

Ar =) Poh bp + > [Arbh btp + Ardo bo] > (4) 
P,O p 


are the isolated (unperturbed) Hamiltonians of the left superconductor, the quantum dot, 
and the right superconductor, respectively. 

We consider only a quantum dot with a negligible intra-dot Coulomb interaction. The 
reason is that we want to focus on the superconducting proximity effect on the QD, which 
is due to Andreev reflections; whereas a large charging effect caused by intra-dot Coulomb 
repulsion suppresses the transport mediated by Andreev process [AB]. Study on the competi- 
tion of a medium intra-dot Coulomb interaction with Andreev reflections in S-QD-S junction 
will be carried on in the coming work. Furthermore, we consider only one single energy level 
inside the QD. The generalization to QDs with several energy levels is straightforward. 

The tunneling term 


Hr = y Wei d + c.c.] + Y [Robi do +c.c.], (5) 


k,o p,o 


describes the electron-transfer between the QD and the leads. 


Current Formula. The current flowing across the junction is [6] 


r=- Tre [7 E lag Wet (w) + GF wet w}, (6) 


where gz, G are the Keldysh Green functions of the left superconductor and the QD in 


Nambu space [A]. The retarded Green function of the superconductor is 


ALR 
Sir (w) = (—im)pryrlw+ 05H |,. e i (7) 
L/R 1 
where pr/R(w) = is the density of states (DOS) of a BCS superconductor. After 
LOS LE 


solving the matrix Dyson equation, we obtain the retarded Green function of the QD: 
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D r r = r 
(w) > (Woy g (w -E (w) 
where 
1 
0 
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is the unperturbed retarded Green function of the QD. The QD’s self-energy 


Zw) =3 81 (u) + Egr (o) , (10) 


describes the proximity effect of superconductors on the QD; here we have used the wide 


Oo band-width approximation [f]. D(w) in Eq. [8] is the determinant of the matrix (g"~'—”): 


= D(w) = [g" OR — =" (w),1][9" (won = Y (wal — 2 (w) X" (0), (11) 


Choice of Integration Contour. The supercurrent formula Eq.|6] can be expressed as 


= Ref [7 Hiag dan 9f + ee (12) 


In order to avoid calculating this integration numerically, we shall transform it into a contour 
integral. First, using the fluctuation-dissipation theorem GS (w) = [G* (w) — @ (w)|f(w), 
where G stands for any Green function in the equilibrium state and f(w) is the Fermi 


distribution function, we have 


he = Ref FIG (00) 9% Wa fle) — Cn Wa fl}. (1 


Then we make a change of the integration variable w > w — i0* for G} (w),. g} (W)o, f(w) 


and w > w+ i0* for G4 (w),, 97 (w)., f(w) to divide Isc into two integrals along different 


paths: 
_4el ee + r “at “+ 
Le Reef gp (Ga (w — i0 Ja gL (w — i0 ),, fw - 00 )] 
BT zu + “yt “yt 
7 er [GS (w + i0+) , Cr (w+ i0 ),, fw + 20 yy. 


Since the above two integrands are actually the same, we define a new integrand equal to 


them: 


J(w) = G3 (w - 10+) „gr (w - 10+), fw - i0*) 
= G4 (w +i0*) „ar (wt i0t), fw + i0*) (15) 
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Thus the current integration can be rewritten into a contour integral in the complex-w plane: 


del +oo-+10+ +oo—40+ E _ 4er 
i= Rell Gi -f ie) re f Hw) (16) 
h oo+i0t oo—i0+ 
The integration path C = f| aa ee Nine =. is a close contour which lies infinitely close to 
the real w axis. The integral can then be evaluated analytically using Cauchy’s theorem. 
This choice of the integration contour has the advantage of leaving all (infinitely many) 


ee arising from f(w) outside the contour, thus avoid the calculation 


2n+1)7 
er), 


the poles at w = 
of the summation over all the Matsubara frequencies (w = i 

Poles and definition of Andreev bound state. Since the current is proportional to the 
sum of the residues of the integrand J(w) inside the contour C, only real poles of J(w) 
contribute. For J(w) = G7 (w — 204), 97 (w — i0), f(w), there are only two pairs of real 
poles: w = +A and w = +e*(|e*| < A). The poles at w = +A originate from the singularities 
of the DOS of the BCS superconductors. 

The poles at w = +e* come from the QD’s Green function and are two real roots of the 
equation D(w) = 0 inside the gap ( |e*| < A ). These two poles give the positions of the 
quasi-particle states of the QD inside the gap. The fact that they are real indicates that 
these two quasi-particle states are exactly bound states, usually called the Andreev bound 
states (ABS). 

Before coupling to the superconductors, the unperturbed QD has a bound state at €o 
for the electron spectrum and one at -—eo for the hole spectrum. After coupled to the 
superconductors, through the electron-transfer with the two superconductors via Andreev 
process, the electron part and the hole part of the QD’s Green function are coupled together, 
and give rise to significant modifications to the QD’s energy spectrum: the original bound 
state at €o for electron (or —ey for hole) renormalizes into two symmetrical bound states at 
+e* for both electron and hole spectrum. 


To express e* as the algebraic function of ey, $ and I’, we transform the equation D(w) = 0 


into a quartic equation of x = € 
2° — (2A? — AT?)z? + (A* + 2A?B? — AT? A’)? + (24?24?Bz + (A?B?)? = 0; 
A =A +e HI’, B= etot Beet, (17) 
This quartic equation of x has only one positive real solution xy less than A?: +e* = +/Zo. 


The quartic equation is easy to solve algebraically; but we omit the lengthy solution. The 


dependence of e* on eo, $ and T is shown graphically in Bj. 


Residues and analytical result of the current. The current is obtained by adding up the 


residues of the four poles: 


ke, 4 
ñas = Im{ Dre sJ(w)l, (18) 


where w; = +e*,+A. The residues of J at e* and —e* are: 


ResSy=+e J (w) = „im [= = = 2a (w — om gi (w _ i0*) 


r 2° f(+e*) sinó 
FR 2 A2 — e*2 2e*[(A? — e2) + (2A? = e*?) + A? sin? e 


Fœ), (19) 
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Im[Re Sw=+e* J(w)] = 


T 
VA 


The residues at +A turn out to be zero since 


1 
DEA] = =0 20 
(+A) DEN (20) 
Thus the supercurrent is simply 
eq? A? sin $ tanh & 
Ls = he A? — &2 l(A2 2 m T2 202 QJ (21) 
e (A? — e?) + Ir (2A O 2 „A sin S| 


where ¢ = ¢@r — $r is the phase difference of the two superconductors. 

Discussion. In addition to the splitting of the original energy level ey into the two 
ABS, another modification to the QD’s energy spectrum is that the QD begins to have 
a continuum energy spectrum outside the gap. In fact, one way to calculate the current 
integral ( I. = SE Rd se [Ga (o) 97 Wa + GF (w)1, 97 ail} ) is to divide the 
QD's spectrum into To continuum part ( |w| > A ) and the discrete part ( |w| < A ), and 
then evaluate their contributions to the current separately. However, this method involves 


a numerical integration and faces the problem of explaining why the peaks at w = +A in 


Icontinuum VS. €o Curve and in Igiscrete VS. €o Curve exactly cancel out when the two parts add 
up to give the total current, see Fig.1. 

In contrast, from our approach, this exact cancellation can be shown naturally. The cur- 
rent is obtained by adding the residues of the four poles. Since the residues at +A vanish, the 
total current curve has no peaks arising from the DOS singularities of the superconductors. 

Now only residues at +e* contribute to the total current and give rise to a sharp peak in 
I vs. €o curve at eg = 0. The reason for the appearance of the peak at ey = 0 is that when 
€o = 0, the upper and the lower ABS have same phase. As e, departs from the Fermi level, 
the phase of the upper and lower ABS begins to differ and this difference severely reduces 
the supercurrent across the junction. Fig.2 shows the central peaks at different I’, the peaks 
broaden when T increases. 

A direct application of this result is a way to align the energy level of the QD to the 
Fermi level of the superconductor in the actual experiment: when changing the gate voltage, 
the state with the maximum supercurrent is at ey = 0. 

Conclusion. In summary, we have calculated analytically the supercurrent across an S- 
QD-S junction. The central results are Eq.[18] and Eq.[21]. They give the analytical result 
of both the positions of ABS and the supercurrent. We show two additional advantages 
of our method: (1) It presents a natural definition of the Andreev bound states inside the 


QD. (2) It explains the disappearance of the peaks at the gap edge in I vs. ey curve. The 


method devised in the paper can also be used in computing other equilibrium properties of 
the superconductor-normal hybrid systems. 

We thank Q.F. Sun for his stimulating remarks. This work is supported by National 
Science Foundation of China under the Grant No. 10074001. 
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FIGURE CAPTIONS 


Fig. 1 Supercurrent vs. eg when T = 0.1 and ¢ = 5. The dash line is from discrete spectrum, 
and the dot line from continuum spectrum. The solid line is total current. Note that 


the continuum spectrum contributes negative current. 


Fig. 2 Supercurrent vs. ey with ¢ = 5 for different I. Different curves are for T = 0.1, 0.3, 


1.0, respectively. Note the broadening effect as I increases. 
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